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We study the effect of diagonal and off-diagonal electron-phonon coupling in the ground state prop- 
erties of the t-J model. Adiabatic and quantum phonons are considered using Lanczos techniques. 
Charge tiles and stripe phases with mobile holes (localized holes) are observed at intermediate 
(large) values of the diagonal electron-phonon coupling. The stripes are stabilized by half-breathing 
modes, while the tiles arise due to the development of extended breathing modes. Off-diagonal terms 
destabilize the charge inhomogeneous structures with mobile holes by renormalizing the diagonal 
coupling but do not produce new phases. Buckling modes are also studied and they seem to induce 
a gradual phase separation between hole rich and hole poor regions. The pairing correlations are 
strongly suppressed when the holes are localized. However, in charge inhomogeneous states with 
mobile holes no dramatic changes, compared with the uniform state, are observed in the pairing 
correlations indicating that D-wave pairing and moderate electron-phonon interactions can coexist. 

PACS numbers: 71.10.Fd, 74.25. Kc, 74.81.-g 



I. INTRODUCTION 

Electron-phonon interactions (EPI) are at the heart 
of the pairing mechanism in the BCS theorji that suc- 
cessfully describes traditional superconductors. It is be- 
lieved, however, that EPI cannot produce the high crit- 
ical temperatures observed in cuprate superconductors 
and, for this reason, most of the research in the field has 
focused on the interaction between charge and magnetic 
degrees of freedom as a source of the pairing mechanism 
in the cupratesi^ 

Despite almost 20 years of intensive research, the pair- 
ing mechanism in the cuprates is still unknown. The 
experimental evidence indicating that phonons and lat- 
tice degrees of freedom play an active role in these 
materials keeps mountingi^iSiSii In addition, inhomoge- 
neous structures described as stripes, patches, tiles, etc., 
have been detected in the underdoped regime of several 
cuprates The emerging phase complexity is reminis- 
cent of the experimental data for manganites where com- 
peting electronic, magnetic and lattice degrees of freedom 
are responsible for the rich phase structure, IS It is also 
likely that the strong electronic correlations present in 
these materials, in particular the proximity to a Mott 
insulating phase in the cuprates, could radically change 
the effects of EPI with respect to those in conventional 
metals and superconductors. 

Most theoretical efforts that tried to incorporate EPI 
in models for the cuprates were performed in the 90 's 
and they focused on whether the interactions could pro- 
duce D-wave pairing (instead of S-wave)ii*i^, and on the 
tendency of phonons towards the stabilization of charge 
density wave (CDW) states that could compete with su- 
perconductivity (SC). In order to observe CDW states 



most of these calculations were performed at quarter 
filling. ^^'M Currently, our aim is to understand whether 
EPI stabilize or destabilize charge stripes, and what kind 
of inhomogeneous textures, if any, develop. Instead of 
focusing on the phonons as an alternative source of the 
pairing mechanism we want to explore whether the in- 
terplay of magnetic and phonon interactions with the 
electrons may lead to an enhancement of the pairing al- 
ready observed in purely electronic models of high Tj, 
superconductors 

In the early studies only diagonal EPI were considered 
and the effects on hopping terms and Heisenberg cou- 
pling due to the lattice distortions (off-diagonal terms) 
were disregarded. Most of the numerical work was per- 
formed with on-site Holstein phonons, and when more 
extended modes were considered the focus was on breath- 
ing modes^'* that would tend to stabilize CDW states and 
in buckling modes that were known to produce the tilting 
of the octahedra in YBCOjii At present, the experimen- 
tal evidence indicates that the electron-phonon coupling 
to the breathing mode is strongly anisotropic and, as a 
result, half-breathing modes seem to have the strongest 
EPI in the high Tc cupratesi^ In addition, some authors 
believe that off-diagonal couplings could be stronger than 
the diagonal ones.^^ 

Most of the recent studies of EPI in models for the 
cuprates have been performed using mean- field, slave- 
boson, or LDA approximationsiiSiiLiSiiS Here, we pro- 
pose to revisit the t-J model and study the effects of EPI 
with unbiased numerical techniques. In order to find the 
relevant phononic modes we will start by studying adi- 
abatic phonons, i.e. uj = 0. Results at finite frequency 
will be presented for the most relevant modes only. 

The paper is organized as follows: Section II describes 
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the study of adiabatic phonons. Diagonal coupling to 
breathing modes is discussed in subsection A while sub- 
section B contains results for diagonal coupling to buck- 
ling modes. The effect of off-diagonal couplings in both 
cases is discussed in subsection C. Quantum phonons for 
the physically important case of half-breathing and buck- 
ling modes are discussed in Section III. Conclusions and 
final remarks are the subject of Section IV. 



all previous exact numerical studies were carried out on 
smaller clusters, -s/S x VS in RefQ and vTO x VTO in 
Ref 13. In addition, larger tilted clusters do not possess 
the symmetry to fit horizontal or vertical stripes. Notice 
that the ionic displacements Wi^^ will be obtained by find- 
ing the configuration that minimizes the total energjiSS 
and without imposing a priori a given pattern for them. 
Details of the procedure were presented in Ref.2Q, and 
will not be repeated here. 



II. ADIABATIC PHONONS 
A. Breathing modes 



The Hamiltonian for the t-J model with adiabatic 
phonons is given by 

H = - ^ tij(cJ^Cj^ +h.C.)+^ JijSi-Sj-fiJe-ph+iJph, 
(ij)'^ (ij> 

(la) 

with 



and 



where 



and 



U,i^t{l + Xt[u(i) + u{i)]}, (3) 



Ji,j = J{l + A,;Ki) + u(j)]}. (4) 



C;^ creates an electron at site i = (ix,iy) of a two di- 
mensional square lattice with spin projection a; Si is the 
spin-i operator at site i, (ij) denotes nearest-neighbor 
lattice sites, t is the hopping amplitude, and J > is 
the antiferromagnetic exchange coupling. In the follow- 
ing, t = 1 and J = 0.4 will be adopted unless other- 
wise stated. Doubly occupancy is not allowed in this 
model. Hi = ^la'^ic is the electronic density on site 
i. Aq is the diagonal electron-phonon coupling and k is 
the stiffness constant. At and Aj are off-diagonal cou- 
pling constants for the hopping and Heisenberg interac- 
tions respectively. The lattice distortions Uj.^ measure 
the displacement along the directions /i = f or y of oxy- 
gen ions located at the center of the lattice's links in the 



equilibrium position, i.e. 



— 0. Ground state prop- 



erties of this model will be studied using the Lanczos 
method^ on square clusters with periodic boundary con- 
ditions (PBC). The results will be presented on 4 x 4 
clusters. Despite the small size we want to remark that 
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FIG. 1: Results for the diagonal breathing mode for four holes 
in a 4 X 4 lattice and J = 0.4. a) Charge structure factor for 
different values of the momentum q as a function of the di- 
agonal electron-phonon coupling. The full (dashed) line indi- 
cates the region where the tile (stripe) ground state is stable; 
b) Pairing correlations. The different symbols represent dif- 
ferent ground states described in the text. 



We will start by discussing the effects of the diagonal 
couphng, i.e., A* = Aj = in Eqs. (3) and (4). For 
the case of four holes, i.e., (n) = 0.75 in a 4 x 4 lat- 
tice, we have observed that the ground state presents 
a uniform charge distribution up to Xq/k = 0.3. The 
probability of hole occupation {nhij)) is the same for ev- 
ery site i and the charge structure factor C(q), shown in 
Fig. da), does not have a sharp peak for any value of the 
momentum q. Notice that the uniform ground state in 
the 4x4 lattice is three-fold degenerate, i.e. states with 
momentum q = (0, 0), (tt, 0) and (0, tt) have the same en- 
ergy. The pairing^ ^ correlations are not the same in these 
three degenerate uniform states. In Fig. ^h) the filled 
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black circles indicate the averaged pairing correlations 
at the maximum distance along the diagonal. The pair- 
ing along the maximum diagonal distance for the state 
with momentum (0, tt) is indicated by the dashed line. 
For Aq/k > 0.3 the uniform ground state is replaced 
by a charge inhomogeneous state with a tiled structure 
schematically displayed in Fig. |2Ia). The size of the cir- 
cles is proportional to the probability of finding a hole in 
the corresponding site. The large circles indicate a prob- 
ability of around 0.3 while the smaller ones correspond 
to about 0.2, which indicates that it is more likely to 
find the four holes in the tiles defined by the big circles 
than outside them. The state is stabilized by "extended 
breathing" modes as denoted by the arrows displayed in 
Fig. Ela), which are proportional to the ionic displace- 
ments Ui^^. Although we are not aware of breathing or 
extended breathing modes being mentioned as relevant in 
the cuprates it is worth mentioning that so called "tiled" 
structures have been observed in STM experiments on 
NaxCa2-xCu02Cl2i^ In addition, all efforts to stabilize 
the observed tiled state in the t-J model without phonons 
have been unfruitful^ The pairing correlations in this 
state, denoted by the triangles in Fig. ^b), are slightly 
suppressed but not destroyed. This state remains stable 
up to Aq/k = 0.62 but from Aq/k > 0.38 it coexists with 
a stripe state schematically shown in Fig.|2Ib). As in the 
previous case, the large (small) circles indicate a site hole 
density of 0.3 (0.2) in the figure. 

In agreement with neutron scattering results for the 
cuprates, four holes in the 4x4 lattice produce two 
stripes as opposed to only one. '' Also, in agreement with 
experimental data^, we observed that vertical (horizon- 
tal) stripes are stabilized by half-breathing modes along 
X (y). The ionic displacements for vertical stripes are 
presented in Fig. |2Ib) . The stripe state has a maximum 
in the charge structure factor at momentum q = (tt, 0) 
(vertical stripes) or (0, tt) (horizontal stripes) which can 
be observed in Fig. ^a). An interesting characteris- 
tic of the stripe state is that the pairing correlations 
at the maximum diagonal distance (filled diamonds in 
Fig. ^b)) slightly increase when compared with the re- 
sult obtained from the average of the three uniform states 
(filled circles). This result is important because it in- 
dicates that D-wave pairing correlations can survive in 
stripe-like states as long as the holes are mobile. 

When the stripe-state is stabilized the triple degener- 
acy in momentum of the ground state is broken. The 
ground state is now a singlet with momentum (0, tt) 
((7r,0)) for vertical (horizontal) stripes. The observed in- 
crease in the pairing correlations does not seem to arise 
due to the break down of the momentum degeneracy since 
the pairing correlations in the uniform state with momen- 
tum (0,7r) (dashed line in Fig. 1(b)) are weaker than in 
the stripe state. We have also observed that in the stripe 
state the pairing correlations are stronger in the direction 
perpendicular to the stripes than in the direction parallel 
to themiSi 

The stripe state becomes the only ground state for 



0.6 < Aq/k < 0.8. The magnetic structure factor S{q) in 
the stripe state has a peak at (7r/2, tt) which indicates a 
TT-shift. 

It is important to remark that a striped ground 
state has been very difficult to stabilize in the t-J 
model without using external fields or special boundary 
conditions. ^'^ However, our results show that a moderate 
diagonal electron-phonon coupling to breathing phononic 
modes stabilizes them in a wide region of parameter 
space. Local lattice distortions in the direction perpen- 
dicular to the stripe, i.e., half-brcathing modes, are re- 
sponsible for this. Coincidentally, these are the modes 
that most strongly couple to electrons in the cupratesi^ 
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FIG. 2: (a) Schematic representation of the tile structure sta- 
bilized by extended breathing phononic modes; the radius of 
the circles is proportional to the hole density Xi = 1 — rn and 
the arrows represent the ionic displacements; (b) Schematic 
representation of the charge stripes stabilized by horizontal 
half-breathing phononic modes; (c) Schematic representation 
of the localized holes that result from large values of the di- 
agonal electron-phonon coupling. 

For Aq/k > 0.8 a state with localized holes, with one 
example schematically shown in Fig. 2(c), becomes the 
ground state. The holes are localized by displacements 
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of the neighboring ions as seen in the figure. The proba- 
bihty of finding a hole in the sites with the large (small) 
circles is almost 1 (0), which indicates that the holes 
are localized. The pairing correlations at long distance 
vanish in this state and they are indicated by the star 
symbols in Fig. f (b). 

For two holes ((n) = 0.875) similar inhomogeneous 
charge structures were observed. As for the case of four 
holes, stripe states, shown in Fig. |3Ja), are stabilized by 
half-breathing phononic modes. As before, the size of the 
circles is proportional to the probability of finding a hole 
in the corresponding site. Two holes induce one stripe 
in the 4x4 lattice in agreement with neutron scattering 
results.^ However, notice the different thickness of the ar- 
rows in Fig.|3Ia) showing that the ionic displacements are 
stronger next to the holes. This state has a complicated 
structure in momentum space and the charge structure 
factors present peaks at (7r/2,0), (7r,0) and (37r/2,0). 
A tile structure with competing energy is also present 
in this case. The corresponding structure is shown in 
Fig. ISI^b) . Notice that although it is more likely to find 
the holes in the tile, the holes are not localized since there 
is a finite probability of finding them outside. The ten- 
dency towards pairing is clear because it is more likely 
that the two holes will be found along the stripe or in- 
side the plaquette. Fully localized holes are observed for 
larger values of the diagonal electron-phonon coupling. 

Qualitative changes are observed when 8 holes are in- 
troduced, i.e. at quarter filling. In agreement with earlier 
studie a^^'^^ we found that the phonon mode that leads to 
the ground state in this case is the full-breathing mode 
that stabilizes a checkerboard charge density wave state, 
as shown in Fig. 4(a) for J — 0.4. The probability of 
finding holes associated to the size of the circles indi- 
cates that the holes are localized. The transition from 
the uniform state to the CDW state as a function of the 
diagonal coupling is characterized by a rapid increase of a 
peak at momentum (tt, tt) in the charge structure factor, 
shown in Fig. 5(a), and for the development of incommen- 
surability, along the diagonal, in the magnetic structure 
factor (see Fig. 5(b)). 

We have observed that for larger values of J, a dimer 
state, shown in Fig. 4(b) for J = 1, develops between 
the uniform and the CDW states. The new phase is 
clearly observed by monitoring the behavior of the charge 
structure factor, shown in Fig. 5(c) for J = f. A peak at 
momentum (7r/2,7r) develops because, as it can be seen 
in Fig. 4(b), the dimer state breaks rotational invariance. 
At the same time a peak at momentum (tt, 0) develops in 
the magnetic structure factor (see Fig. 5(d)). The dimer 
state is stabilized by extended breathing phonon modes 
for intermediate values of the EPI, i.e., 0.5 < Aq/k < 0.9, 
for 3 — \. The dimers are localized pairs. This state was 
not observed in previous studies which considered only 
local breathing modesii2*i^ This result is not surprising 
since it is well known that the t — J model presents phase 
separation at J/t > 3 in the absence of EPF Thus, the 
proximity to phase separation causes the formation of 
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FIG. 3: (a) Schematic representation of the charge stripes 
stabilized by horizontal half-breathing phononic modes for 2 
holes in a 4 X 4 lattice; the radius of the circles is proportional 
to the density of holes and the thickness of the arrows is pro- 
portional to the magnitude of the ionic displacements, (b) 
Same as (a) but for the tile structure stabilized by extended 
breathing phononic modes. 



dimers at moderated values of Xq/k as J increases. For 
a fixed J the checkerboard CDW state is stabilized at 
sufficiently large values of the EPF 

The pairing correlations in the CDW and dimer states 
vanish at large distances which is not unexpected since 
in both cases the holes are localized. 

In summary, the previous results show that the diag- 
onal electron-phonon coupling to breathing modes sta- 
bilizes charge inhomogeneous structures, some of which, 
like stripes and tiles, have been observed in the cuprates. 
The pairing correlations are depressed when the holes are 
localized. However, in charge inhomogeneous states with 
mobile holes the D-wave pairing correlations survive. Al- 
though it is not possible to infer that long range correla- 
tions will develop in larger lattices, these results indicate 
that the interaction between magnetic and phononic de- 
grees of freedom does not necessarily destroy pairing. 



B. Buckling modes 

For completeness, the effect of buckling modes will also 
be considered. Buckling modes are believed to be impor- 




FIG. 4: (a) Schematic representation of the CDW state sta- 
bihzed by breathing phononic modes using 8 holes in a 4 x 4 
lattice for J = 0.4; (b) Schematic representation of the dimer 
state stabilized by extended breathing phononic modes at 
quarter filling (8 holes) for J = 1. 




FIG. 5: (a) Charge structure factor for 8 holes in a 4 x 4 
lattice for J = 0.4; (b) Magnetic structure factor for the same 
parameters as (a); (c) Same as (a) for J = 1; (d) same as (b) 
for J = 1. 




tant in YBCOi242& and it was suggested that in combina- 
tion with the magnetic interactions they could produce 
D-wave electronic pairing jii The buckling mode is incor- 
porated by replacing u{i) in Eq.(l) by 

-"(i) = Ui^x + Ui-x.j; + + Ui-y,a, (5) 

while ^i.j and Jij have to be replaced by 

ti,\+^^ = t[l + XtUi^^], (6) 

and 

Ji,i+^ = J[l + AjMi_^]. (7) 

The most remarkable effects of adiabatic buckling 
modes are the following: 1) stripe and tile states quali- 
tatively similar to the ones observed with half-breathing 
modes are stabilized in a very narrow region at interme- 
diate values of the diagonal coupling, and 2) the ground 
state phase separates for large values of Aq/k. 

We have observed the phase separated state by dop- 
ing with 2, 4 and 8 holes. In the case of 2 holes pair 
localization appears. However, for 4 holes it can be seen 
that the holes localize in clusters instead of remaining 
isolated as it was the case with breathing modes (see 
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FIG. 6: Schematic representation of the phase separated 
states induced by buckling phononic modes for 4 holes in a 
4x4 lattice: a) charge stripes; (b) the tile structure. 



Fig. 2(c)). Examples of the phase separated clusters for 
4 holes arc shown in Fig. 6. The stripe (tile) structure 
in Fig. 6(a) (Fig. 6(b)) resembles the stripe (tiles) in 
Fig. 3(a) (Fig. 3(b)) obtained for 2 holes with breath- 
ing modes. However, we can see that the probability of 
hole occupation is almost 1 in the stripe and the tile in- 
homogeneities of Fig. 6 indicating localization. In the 
phase-separated state the holes are localized by strong 



uniform displacements of the oxygens. In the case of the 
stripe, only the oxygens in the links along the stripe are 
strongly displaced while the tile is stabilized by a strong 
displacement of the oxygens in the links connecting the 
holes. In the interval 0.75 < Aq/k < 2.5 the tile and the 
stripe states have very similar energies. 




FIG. 7: (a) Charge structure factor as a function of the di- 
agonal electron-phonon coupling for the buckling mode and 2 
holes in a 4 X 4 lattice; (b) Pairing correlations at the maxi- 
mum distance r-max for the same parameters as (a) 



At quarter filling, i.e., 8 holes, it can be seen that the 
holes form a cluster so that two neighboring rows (or 
columns) in the 4x4 cluster have electrons and the other 
two rows (or columns) only have holes. Again, the holes 
are localized by strong displacements of the surrounding 
oxygens. As expected, the pairing correlations immedi- 
ately vanish in the phase separated state. 



Between the uniform and the phase separated states, a 
region with charge inhomogeneous structures is observed. 
For two holes we observe in Fig. 7(a) that a peak develops 
in the charge structure factor for Aq/k — 0.25. Tiles 
and stripes similar to those displayed in Fig. 3 become 
the ground state. As the couphng increases the holes 
become more and more localized in the stripe or the tile 




FIG. 8: (a) Charge structure factor as a function of the di- 
agonal electron-phonon coupling for the buckling mode and 4 
holes in a 4 X 4 lattice; (b) Pairing correlations at the maxi- 
mum distance Vmax for the same parameters as (a). 



but phase separation occurs suddenly for Aq/k > 1 when 
the holes form a localized pair. In Fig. 7(b) it can be seen 
how the pairing correlations are reduced in the charge 
inhomogeneous phases as localization increases. 

A similar behavior, displayed in Fig. 8, is observed for 
4 holes. The remarkable result is that buckling modes 
appear to be able to stabilize the same kind of stripes 
than the half-breathing modes produce, but only in a 
very narrow region. Let us study in detail what is the 
underlying ionic configuration. In Fig. 9(a) we observe a 
schematic representation of the stripe state for Aq/k = 
0.4. The hole density distribution is similar to the one 
for half-breathing modes shown in Fig. 2(b) but the ionic 
displacements are not the same. Clearly there is a break 
down of the rotational invariance so we could call the 
mode "half-buckling" . We notice that only the ions in the 
direction parallel to the stripes are displaced, while the 
breathing modes stabilize the stripe via displacements 
of the ions in the links perpendicular to the stripe. The 
oxygens in the links that connect the hole-rich sites move 
in one direction while those in the links that connect the 
hole poor sites move in the opposite direction. 

This phononic configuration appears as an alterna- 
tive to the half-breathing modes for stabilizing dynamic 
stripes but it is very unstable. The induced charge inho- 
mogeneity is very small, of the order of 1%, and as we 
will see below, the state is rapidly replaced by a precursor 
of the phase-separated state. 

For Aq/k = 0.4 a tile state is also stabilized. It is 
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shown in Fig. 9(b). Wc observe that only the oxygens 
in the links that connect the sites that form each tile 
are displaced. The oxygen displacements are opposite in 
hole rich and hole poor regions. It seems as if a depres- 
sion in the position of the oxygen ions tends to localize 
the holes, like a depression in a mattress would tend to 
localize marbles. It is interesting to monitor how the 
stripe and tile configurations evolve to the phase sepa- 
rated states shown in Fig. 6. At Aq/k = 0.5 the stripes 
become distorted. An upward displacement of the oxy- 
gen ions along the last column in Fig. 10(a) effectively 
reduces the hole density in that column. A similar ef- 
fect is observed in the tile state where, as we can see in 
Fig. 10(b), an upward displacement of the oxygens in one 
of the tiles pushes the holes out. We have observed that 
for 0.4 < Aq/k < 0.5 the tile state has energy slightly 
lower than the stripe state. 




FIG. 9: (a) Schematic representation of the distorted charge 
stripes with weakly localized holes stabilized by buckling 
phononic modes for 4 holes in a 4 x 4 lattice at A§/k = 0.4; 
(b) Same as (a) but for the tile structure. 



An important qualitative change occurs at \1/k = 0.6. 
Both in the stripe (Fig. 11(a)) and in the tile (Fig. 11(b)) 




FIG. 10: (a) Schematic representation of the charge stripes 
stabilized by buckling phononic modes for 4 holes in a 4 x 4 
lattice; (b) Same as (a) but for the tile structure stabilized by 
buckling phononic modes. 



configurations the holes get localized in the region of the 
lattice in which the oxygens are depressed. The stripe 
state has the lowest energy in this case. 



Finally, at Aq/k = 0.7 a structure similar to the one 
shown in Fig. 6 develops. The ion displacements become 
stronger as the diagonal coupling increases as well as the 
local differences in hole density. Actual phase separa- 
tion, with n{i) « or 1 in every site, i.e., with the holes 
completely localized, appears to develop for Aq/k > 1. 

In Fig. 8(b) we can see the behavior of the pairing 
correlations in the different phases. It is interesting to 
compare with the corresponding results in Fig. 1(b) for 
the half-breathing case. The pairing correlations in the 
uniform state are, of course, the same as in Fig. 1(b). 
Thus, to simplify the figures we only present results for 
the state with momentum (0,7r) in the imiform case de- 
spite the fact that we know that the ground state is triply 
degenerate in momentum. The filled circles indicate the 
pairing along the maximum diagonal distance. 

When the double stripe is stabilized for Aq/k = 0.4 





FIG. 11: (a) Schematic representation of the charge stripes, 
precursor of the phase separated state, stabilized by buckling 
phononic modes for 4 holes in a 4 x 4 lattice for Xq/k, = 0.6; 
(b) Same as (a) but for the tile structure. 



we observe a behavior very similar to the one for stripes 
in the breathing case. However, the stripes are destabi- 
hzed by a phase which is a precursor of the phase sepa- 
rated state where the holes are localized by the buckling 
modes (see Fig. 10). The area of the localization region 
decreases with the diagonal EPI until the holes form a 
localized cluster and phase separation is reached. The 
pairing correlations, as expected, decrease sharply. 

An interesting detail is that, as in the half- breathing 
case, in all stripe states the pairing correlation is larger 
in the direction perpendicular to the stripes that in the 
direction parallel to them. 

At quarter filling the transition from the uniform to 
the phase-separated state is more abrupt. In fact, only 
for A§/k = 0.9 the localization region is larger than two 
rows or columns. In Fig. 12(a) it is shown how the weak 
peak at (tt, tt) in the charge structure factor characteris- 
tic of the uniform ground state is replaced by a strong 
peak at (0,7r/2) indicating that half the lattice is filled 
with holes and the other half is empty. On the other 
hand, the magnetic incommensurability characteristic of 
the uniform quarter-filled state dissapears in the phase 
separated region as it can be seen in Fig. 12(b). 
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FIG. 12: (a) Charge structure factor as a function of the 
diagonal electron-phonon coupling for the buckling mode and 
8 holes in a 4 X 4 lattice; (b) Magnetic structure factor for the 
same parameters as (a). 



Summarizing, the buckling modes induce a lattice de- 
pression where the holes concentrate. A "half-buckling" 
mode that plays a role similar to the "half-breathing" 
mode is rapidly destabilized by the diagonal coupling in 
favor of a state with localized clustered holes and sup- 
pressed pairing correlations, which is a precursor of the 
phase-separated state. 



C. Off-Diagonal Electron-Phonon coupling 

At the adiabatic level, we have observed that the ef- 
fect of the off-diagonal terms is to destabilize the charge 
inhomogeneous states with mobile holes that appear in 
between the uniform and the localized holes or phase sep- 
arated phases. This can be seen in Fig. 13 where a phase 
diagram in the plane At vs Aq/k for breathing modes is 
presented for 4 holes and Aj = A^. It can be seen that for 
Aq/k w 0.7 the charge inhomogeneous stripe/tile state is 
replaced by an uniform ground state for At = 0.1. As 
Ao increases, larger values of At are needed to restore 
the uniform ground state. The region labeled M in the 
figure is characterized by charge inhomogeneous states 
with mobile holes. This phase is overcome by the uni- 
form state for larger values of the off-diagonal coupling. 
We have not observed any enhancement of the pairing 
correlations due to the off-diagonal coupling. 

For buckling modes, see Fig. 14, we observe that small 
values of At just renormalize the diagonal coupling to a 
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smaller value so that the uniform ground state remains 
stable up to larger values of Aq/k. However, larger values 
of Xt destabilize the inhomogeneous phase with mobile 
holes labeled Ad in the figure, and for A* > 0.2 the sys- 
tem goes from a uniform to a phase separated state as a 
function of Aq/k. 

The off-diagonal coupling does not generate new 
phases, it simply seems to rcnormalize the diagonal cou- 
pling. In fact, its effect is to destabilize the inhomoge- 
neous phase with mobile holes identified by M in Figs. 13 
and 14. This phase is more easily destabilized for buck- 
ling than for breathing modes. 
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FIG. 13: Phase diagram for 4 holes as a function of the diag- 
onal EPI Aq/k and the off-diagonal coupling At for breathing 
modes. In this case \j = At. M characterizes charge inhomo- 
geneous states with mobile holes described in the text. 



III. QUANTUM PHONONS 

The results presented in the previous section indicate 
that the phonon mode that most closely reproduces the 
experimental results observed in the cuprates is the half- 
breathing. In addition, neutron scattering experiments 
show that this is the mode the most strongly coupled to 
the charge in the high Tc cuprates. For this reason we 
decided to study the effects of the half-breathing mode 
beyond the adiabatic approximation. 

In order to study the electron-phonon interaction at 
finite frequency the electron-phonon part of the Hamil- 




FIG. 14: Phase diagram for 4 holes as a function of the diag- 
onal EPI Aq/k and the off-diagonal coupling At for buckling 
modes. In this case Aj = At. M characterizes charge inhomo- 
geneous states with mobile holes described in the text. 



Ionian (Eq.(lb)) should be replaced b}ii^ 

where the minus sign between the densities corresponds 
to the breathing mode while the plus sign corresponds 
to the buckhng mode that will be discussed later. b\ ^ 
creates a phonon at the link i, /i with frequency lo and 

Aq = g\J 9;^, where g is the electron-phonon coupling 
and ni is the ionic mass. Eq.(lc) should be replaced by 

H,u^u;Y^{b\j,^, + \). (86) 

The off-diagonal effects are given by 

ti,j=ni + At[i?(i)+B(j)]}, (9) 

and 

Jij = J[l + Aj[i3(i) + B(j)]}, (10) 

where 

(11) 

In order to simplify these expressions we make the Fourier 
transform of phonon and electron operators. Then, in 
momentum space, the interaction term becomes: 

k.q,/i,cj 

(12) 

The effective electron-phonon coupling constant 
A(q, ^) = Ao(l T e~''i '^) is clearly anisotropic. For 
breathing modes (negative sign inside the parenthesis) 
it is strongest for cimax = (t^, tt) and, thus, in the 
simplest approximation only the phonon mode 6q„„^ is 



maintained li^ii^ For half-breathing modes qmax = (tt, 0) 
or (0, tt). Under this approximation the Hamiltonian for 
half-breathing modes along x becomes: 



rrbr 



2Ao(6t+6)^ni(-ir 



where Xi is the x-coordinate of site i. 

1, 



rrbr 



Lu{b^b+-). 



(13) 



(14) 



The ofF-diagonal effects on the hopping and the Heisen- 
berg coupling vanish along the direction of the half- 
breathing modes (i in this case). Along the perpendicu- 
lar direction they are given by 



t 



t[i+4Xt{-ir'ib^ + b)], 



and 



^U+y-^[l + 4Aj(-l)-'(6t+6)]. 
This behavior was also remarked in RefEE 



(15) 
(16) 




FIG. 15: (a) Magnetic structure factor for different values of 
the momentum as a function of the diagonal electron-phonon 
coupling, for quantum half-breathing modes in a 4 x 4 lattice 
for 2 holes, (b) Charge structure factor, pairing at maxi- 
mum distance, and average number of phonons as a function 
of the diagonal electron-phonon coupling, for quantum half- 
breathing modes in a 4 x 4 lattice with 2 holes. 
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FIG. 16: (a) Magnetic structure factor for different values of 
the momentum as a function of the diagonal electron-phonon 
coupling, for quantum half-breathing modes in a 4 x 4 lattice 
for 4 holes, (b) Charge structure factor, pairing at maxi- 
mum distance and average number of phonons as a function 
of the diagonal electron-phonon coupling, for quantum half- 
breathing modes in a 4 x 4 lattice for 4 holes. 




FIG. 17: (a) Magnetic structure factor for different values of 
the momentum as a function of the diagonal electron-phonon 
coupling for quantum half-breathing modes in a 4 x 4 lattice 
for 8 holes, (b) Charge structure factor, pairing at maxi- 
mum distance and average number of phonons as a function 
of the diagonal electron-phonon coupling for quantum half- 
breathing modes in a 4 x 4 lattice for 8 holes. 
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As in the previous sections, J will be fixed to 0.4; lo 
will be set to 1. 

Let us discuss first the effects of the diagonal coupling, 
i.e.. At = A,/ — 0. An important difference between the 
adiabatic and the quantum case is that with adiabatic 
phonons the ionic displacements were unconstrained and, 
as a result, we observed that in some regions of pa- 
rameter space extended breathing modes with rotational 
symmetry that stabilized tile, CDW, and localized holes 
structures, were energetically more favorable than half- 
breathing modes. With quantum phonons, we will study 
half-breathing modes only and, thus, all the charge inho- 
mogeneous states that will arise will be stripe-like. 

In Fig. 15 results for 2 holes doped in a 4 x 4 lattice are 
presented. In panel (a) we observe the magnetic struc- 
ture factor for different values of the momentum. It can 
be seen that the value at (tt, tt) is always maximum and, 
as a result, magnetic incommensurability does not de- 
velop. The charge structure factor (open circles in panel 
(b)) indicates that a stripe- like structure starts to develop 
at 2Ao = 0.5, i.e., when the average number of phonons 
(star symbols) becomes different from zero. However, 
the pairing correlations at the maximum diagonal dis- 
tance (indicated by triangles) always decrease for finite 
EPI and actually vanish in the stripe phase. This oc- 
curs because the stripe pattern that results in this case 
is different from the physically "correct" stripe shown in 
Fig. 3(b) that we obtained in the adiabatic case. In order 
to obtain this pattern with quantum phonons, finite val- 
ues of (7(7r/2,0) should be consider in addition to finite 
values of 5(7r,0), as explained in Section II. The single 
coupling at momentum (tt, 0) stabilizes a double stripe 
pattern with one hole per stripe. The charge difference 
between the background and stripes is large enough to 
develop a peak at (7r,0), but the behavior of the mag- 
netic structure factor indicates that it does not produce 
a TT-shift. 

In Fig. 16, results for different values of Aq on a 4x4 lat- 
tice doped with four holes are presented. In this case we 
obtain results qualitatively very similar to those in the 
adiabatic case (Fig. 1). The magnetic structure factor 
displayed in Fig. 16(a) develops a peak at q = (7r/2,7r) 
for 0.2 < 2Ao < 0.5, which indicates that a state with 
magnetic incommensurability is stabilized. The charge 
structure factor shown in Fig. 16(b) (with open circles) 
indicates that charge order consistent with two vertical 
stripes develops in the magnetically incommensurate re- 
gion. The uniform ground state is replaced by an inho- 
mogeneous one when the number of phonons (star sym- 
bols) becomes finite at 2Ao > 0.2. This is similar to the 
effect observed with adiabatic phonons. For 2Ao > 0.5 
the charge structure factor results indicates that a static 
stripe-like structure has developed. A 7r-shift does not 
occur since S{(\) has a maximum at (tt, tt). 

A small enhancement in the pairing correlations at 
maximum diagonal distance (triangles in Fig. 16(b)) oc- 
curs as a function of Aq in the uniform phase. As in the 
adiabatic case, D-wave pairing correlations are not dra- 



matically affected in charge inhomogeneous states with 
mobile holes. 

Fig. 17 shows the results at quarter-filling. Instead of 
a CDW stabilized by full-breathing modes, which is what 
we observed in the adiabatic case, the half-breathing 
modes still generate a two-stripe structure. This oc- 
curs because we are still studying the effects of half- 
breathing modes, while phonon operators with momen- 
tum (tt, tt) would induce the CDW state. The strong 
peak in the charge structure factor C at (tt, 0) (open 
circles in Fig. 17(b)) underlines the large charge inho- 
mogeneity which induces magnetic incommensurability, 
as it can be seen from the magnetic structure factor's 
behavior displayed in Fig. 17(a). 
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FIG. 18: (a) Charge (C) and Magnetic (S) structure factors 
for different values of the momentum, as a function of the 
off-diagonal electron-phonon couphng At for quantum half- 
breathing modes in a 4 x 4 lattice for 4 holes. Aj = 2At and 
the diagonal coupling is given by 2Ao = 0.2. (b) Pairing at 
maximum diagonal distance, and average number of phonons 
as a function of the off-diagonal electron-phonon coupling for 
quantum half-breathing modes for the same parameters as in 
part (a). 



The tendency to incommensurate order increases as 
the doping is increased from 2 to 8 holes, as it can be 
seen in Figs. 15, 16 y 17 ((a) panels). This is under- 
standable since the antiferromagnetic ordering between 



the two undoped domains outside the stripes is reduced 
when the stripes become increasingly populated by holes. 

Let us now consider the effects of off-diagonal electron- 
phonon coupling. We will start with the case of 4 holes 
at 2Ao = 0.2 which, according to Fig. 16(b), provides the 
state with the strongest pairing correlations. We have 
set A J = 2Xt^ The state at At = is uniform and, 
like in the adiabatic case (Fig. 1), is triply degenerate in 
momentum. The off-diagonal EPI breaks the degener- 
acy in momentum and selects de ground state with mo- 
mentum (0, tt). As Xt increases the peak in the charge 
structure factor at momentum (tt, 0) remains almost un- 
changed (circles in Fig. 18(a)). The peaks at momentum 
(tt, tt) (squares) and (7r/2,7r) (diamonds) in the magnetic 
structure factor remains essentially unchanged as well. 
While this happens the pairing correlations and the av- 
erage number of phonons decrease shghtly (Fig. 18(b)). 
This trend continues up to AXt = 0.2. Beyond this value, 
there are no indications of the formation of a charge in- 
homogeneity but there are some changes in the magnetic 
channel with an incipient order with momentum (tt, 7r/2). 

We will also study the effects of off-diagonal interac- 
tions for 2Ao = 0.4, i.e., in the region with dynamical 
stripes. The effect of A* in this case is to select the state 
with momentum (0,0) as the ground state. At At = the 
ground state presents a charge inhomogeneous state with 
two stripes according to the large peak at (tt, 0) in C(q) 
(circles in Fig. 19(a)) exhibiting magnetic incommensu- 
rability. The pairing is very reduced. As At increases 
up to 4At =0.1 the charge inhomogeneity is reduced 
and the pairing correlations increase indicating that At 
is renormalizing Aq making it effectively smaller and sta- 
bilizing the charge uniform state. This is clearly seen at 
4At = 0.15 when the ground state changes momentum 
to (tt, 0) and we observed a ground state with properties 
similar, including the pairing, to the case Aq = 0.2 and 
At = 0. The off-diagonal coupling just seems to make the 
diagonal one weaker. 



This level crossing is related to the disappearance 
of the magnetic incommensurability at (7r/2,7r). Once 
this incommensurability has disappeared, a second level 
crossing takes place at 4At ~ 0.225 and the ground state 
momentum is again (0,0). By further increasing At, the 
charge ordering is reestablished. This complex behavior, 
together with the one previously found for 2Ao — 0.2, 
illustrates the interplay between the diagonal and off- 
diagonal EPI. In particular, due to a nonzero Aj, the spin 
degrees of freedom are now directly coupled to the ion 
displacements, and hence magnetic orderings are more 
relevant to the overall behavior than with purely diago- 
nal EPI. 

For completeness we are going to comment on the ef- 
fects of the diagonal electron-phonon coupling to buck- 
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FIG. 19: (a) Charge and magnetic structure factor for 
different values of the momentum as a function of the 
off-diagonal electron-phonon coupling At for quantum half- 
breathing modes in a 4 x 4 lattice for 4 holes. Aj — 2Xt and 
the diagonal coupling is given by 2Ao — 0.4. (b) Pairing at 
maximum diagonal distance, and average number of phonons 
as a function of the off-diagonal electron-phonon coupling for 
quantum half-breathing modes for the same parameters as in 
part (a). 



ling modes with quantum phonons. 

Our studies with adiabatic phonons showed that in 
a very narrow region of parameter space stripes were 
stabilized by a buckling mode that breaks rotational 
invariance. The mode that stabilizes vertical stripes 
would have a maximum coupling constant at momentum 
(tt, 0). Selecting this momentum in Eq.(12) (with the plus 
sign) we arrive to the same Hamiltonian as in the half- 
breathing case, except that the single phonon modes are 
now along the y direction (parallel to the stripe) rather 
that in the direction x as in the half-breathing case, in 
agreement with the results obtained in the adiabatic case 
(Fig. 9(a)). In addition, off-diagonal effects vanish in this 
case. Thus, the same physical properties as in the diag- 
onal half-breathing case are obtained. 

In Ref. it was shown at the mean-field level that a 
diagonal coupling was sufficient to enhance D-wave pair- 
ing correlations. Selecting the plus sign in Eq.(12) we 
observe that the effective coupling constant is maximum 
for q = (0, 0) and vanishes at momentum (tt, tt) in agree- 
ment with the result of Nazarenko et alii. Regretfully 
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keeping only momentum q = (0, 0) in Eq.(12) leads to 
Hl^^^2XoY,{bl + b^)Y,n,, (17) 

i 

and 

H'^i, = LuJ2iblb, + l), (18) 

which does not couple the electrons to the phonons since 
the number of electrons is fixed in our simulations and 
J2i is constant for all states. A study of this case con- 
serving more than two phonon modes would be desirable 
since mean field calculationa^i and results for 2 holes dop- 
ing obtained with a truncated approach^* provided en- 
couraging results towards D-wave pairing stabilization. 

IV. CONCLUSIONS 

Summarizing, we have studied the efi^ects of electron- 
phonon interactions in the t-J model for a variety of 
phononic modes, couplings, and electronic densities us- 
ing small clusters solved exactly. The most remarkable 
result is that charge inhomogeneous states with mobile 
holes are stabilized by these interactions. In particular, 
stripe states which are very difficult to stabilize in the 
t-J model are clearly observed. Tile structures, detected 
experimentally in some cuprates, are stabilized for the 



first time in the t-J model. Our calculations also confirm 
that the half-breathing mode that stabilizes the stripes 
is the most energetically favorable in some regions of pa- 
rameter space. We have also found that buckling modes 
can generate stripe and tile states in a very narrow range 
of electron-phonon coupling constants. Thus, this mode 
also has to be considered in the interpretation of ex- 
perimental data. The pairing correlations vanish in the 
regime of localized holes induced by the EPI but it is 
encouraging that not drastic effects occur in the states 
with mobile holes. In fact, small enhancements have 
been detected in some cases. This should be contrasted 
with the pair breaking effects of the full breathing mode 
in D-wave superconductors. As previously observed in 
the context of the spin-fermion model^** and diagramatic 
calculations f2i diagonal EPI induce charge localization 
while off-diagonal EPI encourages hole mobility. 

Our results indicate that electron-phonon interactions 
need to be considered in order to understand the proper- 
ties of the high-Tc cuprates, in particular, the existence 
of charge inhomogeneous ground states. 
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